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ABSTRACT 

We develop a simple analytical model that tracks galactic metallicities governed by star formation 
and feedback to gain insight from the observed galaxy stellar mass-metallicity relations over a large 
range of stellar masses and redshifts. The model reveals the following implications of star formation 
and feedback processes in galaxy formation. First, the observed metallicity relations provide a strin¬ 
gent upper limit for the averaged outflow mass-loading factors of local galaxies, which is ~ 20 for 
M* ~ 10 ®Mq galaxies and monotonically decreases to ^ 1 for M* ^ lO^^Mg galaxies. Second, the 
inferred upper-limit for the outflow mass-loading factor sensitively depends on whether the outflow 
is metal-enriched with respect to the ISM metallicity. If half of the metals ejected from SNe leave 
the galaxy in metal-enriched winds, the outflow mass-loading factor for galaxies at any mass can 
barely be higher than ~ 10, which puts strong constraints on galaxy formation models. Third, the 
relatively lower stellar-phase to gas-phase metallicity ratio for lower-mass galaxies indicate that low- 
mass galaxies are still rapidly enriching their metallicities in recent times, while high-mass galaxies 
are more settled, which seems to show a downsizing effect in the metallicity evolution of galaxies. 

The analysis presented in the paper demonstrates the importance of accurate measurements of galaxy 
metallicities and the cold gas fraction of galaxies at different redshifts for constraining star formation 
and feedback processes, and demonstrates the power of these relations in constraining the physics of 
galaxy formation. 


1. INTRODUCTION 

Star formation and galactic outflows triggered by feed¬ 
back from star formation are considered the most im¬ 
portant process es in galaxy formation (|Bensonl l201ft 
iMo et all 1201011 . Galaxy surveys find that galaxies 
comprise only a small fraction of baryonic matter in 
the U niverse le.g. iBrinchman n fc Ellisll2000t iCole et ^ 
" ^ a].l 1200311 ■ 


120011 : iDickinson et al.l I200.?I iBell et al.l 120031 1. Assum- 
ing galaxies form in Co ld Dark Matter (CDM) halos 
(jBlumenthal et ahl 119841 , the baryon mass fraction in 
halos hosting a galaxy smaller than the Milky Way 
must decrease rapidly for decreasing halo mass in or¬ 
der to explain the shallow slope of the low-mass end of 
galaxy mass functi o n, pointed o u t in pioneering works b y 
iFrenk et~ni (|1988ll : [Cold (jlQOlh : [White fc Frenkl (1199111 


and elaborated by recent analysis (e.g. iPanastergis et M 


ml 


Behroozi et al.l 12511 iMoster et al.ll2013t iLu et ~ 


2015bll . Strong outflows have been proposed as the 


most important proces s responsible for keeping low-mass 
galax ies baryon poor (jPekel fc SilkI 119861 : iLacev fc SilkI 
EMI e.g .). Many galaxy formation models success¬ 
fully reproduce the number density of low-mass galax¬ 
ies by invoking strong outflows (e.g. [Benson et al.ll2003l : 
[Somerville et al.[[200^ ICuo et al.[[2011[l . Using model in - 
ference techniques, [Lu et al.[ ([20141 and [B enso i (Iml 
found that to match the data the outflow rate needs to 
be more than 10 times higher than the star formation rate 
for low-mass galaxies. Moreover, these outflows must be 
sustained, or at least be recurrent, over cosmological time 
scales. Although strong feedback in starburst galaxies is 
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often observed, it is not clear if the duty cycle of the out¬ 
flow is sufficiently high to be consistent with what seems 
to be required in the models. In spite of the success of the 
assumption of strong outflows in reproducing many sta¬ 
tistical properties of the galaxy population, observational 
evidence for such strong outflows in local galaxies is sti ll 
lacking (e.g. [Bouche et al.l [201^ [Kacnrzak et al.l 120141 . 
Therefore, it is important to seek other independent ob¬ 
servational tests to constrain the strength of outflows. 

Metals other than the primordial species can only be 
produced by star formation. Therefore, metallicities of 
the baryonic matter in different phases in a galaxy are 
expected to put interesting constraints on star formation 
and outflow. For a given galaxy, the total stellar mass 
provides the total metal mass budget to be distributed 
anywhere associated with the galaxy (including in out¬ 
flows). The total amount of metal mass remaining in a 
galaxy, therefore, contains critical information about how 
outflows work. Moreover, how metals are partitioned be¬ 
tween different phases of baryonic matter may give some 
insight into the nature of star formation. Using the com¬ 
bination of gas- and stellar-phase metallicities, one can 
attempt to infer the broad nature of the star formation 
history in a galaxy. 

To extract these pieces of information from observa¬ 
tional data, we develop an analytic model that can logi¬ 
cally connect star formation and outflow processes with 
the observed metallicity relations. Using this model, we 
attempt to gain insight into these metallicity relations 
and to constrain the strength of outflow and the metal¬ 
licity of past star formation. We find that with simple 
but plausible assumptions such a model can match the 
observed metallicity relations and draw inferences about 
star formation and outflow based on existing data. 

The paper is organized as follows. In §|2l we describe 
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the model of the evolution of metal content of galax¬ 
ies which is based on minimal assumptions. We develop 
two different approaches to constrain the mass-loading 
factors of galactic outflows and star formation in differ¬ 
ent environments using data on the mass-metallicity rela¬ 
tions. In §[3l we describe the observational data adopted 
in this paper. We then demonstrate in §2] how obser¬ 
vational data on the stellar- and gas-phase metallicities 
constrain the outflow mass-loading factor and star for¬ 
mation. In §[5l we discuss the conclusions and implica¬ 
tions of these results. Throughout the paper, we assume 
Solar metallicity Zq = 0.0134 (jAsphind et al.l I2009H 
and Solar oxygen abundance 12-1- log(0/if)fn = 8.69 
(|Allende Prieto et al.l 120011: lAsnlund et al.l 1200911 . We 
adopt a Chabrier IMF ( Chabried 1200311 f^ all data and 
modeling. We have converted all the adopted data based 
on these assumptions. Another important quantity we 
need to hx is the chemical yield, y, which is defined as 
the ratio between the mass of newly produced metals 
that are ejected into the interstellar gas and the mass 
locked in long-lived star s for a single stell a r pop ulation, 
as originally defined bv ISearle fc Sargeiitl (|1972ll . Note 
that this is different from another widely adopted defini¬ 
tion of yield in the literature, p, which is defined as the 
mass of newly produced metals per unit gas mass that 
is turned into stars (i.e. these two definitions differ due 
to the recycling of gas from short-lived stars). In this 
paper, we choose to use a rather high metal yield, y = 
0.07, which corresponds to p = 0.038 for the Chabrier 
IMF. This value is at the higher end of values pre¬ 
dicte d by nucleo-synthesis m o dels (iWooslev fc Weaver] 
II995I: IWooslev fc Heeerl 120021 : iNomoto et al.l Il997bllar ). 
We choose to use this high yield to have a conservative 
estimate of the upper limit of the inferred outflow mass¬ 
loading factor. 


2. AN ANALYTIC MODEL FOR GALAXY METALLICITY 
The origin and implication of galaxy metallicity 


Matteucci fc Tornambe 


Peeples fc Somervilld 


fe.ff. iTin 
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119781 

19871: 
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lErbl 120081: IDave et al. 

20121: 
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iiilil: 


Zahid et al.l 12014), sem i -analy t ic galaxy formation 


models (e.g. ~le et al.l l2000t iDe Lucia et al.l 12004 


iNagashima et al.1 120051 : lYates et ahT 2013fl. and hydro- 
dynamical simu l ations (e.g. iFinlator fc Davd 120081: 
Dave et ^120111: iKobavashi et al.ll2007t IWiersma et al.1 


20091 : iMa et al.l 201511 . In this paper, we take the ana¬ 


lytical approach to build a simple model that connects 
the buildup of metal mass in a galaxy and outflows of 
material from that same galaxy. In this model, star 
formation happens in the interstellar medium (ISM), 
which is continuously enriched with metals by star 
formation. We assume that metals are instantaneously 
produced by star formation, and the time delays between 
the nuclear synthesis of different heavy elements are 
ignored for simplicity (i.e. the usual “instantaneous 
recycling approximation”). In addition, we assume that 
the mixing of metal mass in the ISM is perfect. For 
simplicity, we further assume that the newly accreted 
gas does not contribute any signihcant metal mass. 

Under these assumptions, we can write down a set 
of equations to follow the metal mass. The change of 


the metal mass in the ISM when dM* of stellar mass is 
formed can be written as 

dMz_g = j/dM* — ZgdM* — ^ , (I) 

where denotes the metal mass in cold gas, M* 

denotes the mass of long-lived stars, y is the chemical 
yield, dehned as the ratio between the mass of newly 
produced metals and th e mass locked in long-lived stars 
(|Searle fc Sargentlll972h . Zg is the gas-phase metallicity 
at the time of star formation, R is the mass fraction that 
is returned into ISM from short lived stars and stellar 
wind, and r] is the galactic outflow mass-loading factor, 
which is defined as the ratio between the outflow mass 
flux and the instantaneous star formation rate. The first 
term in the right side of the equation represents the metal 
mass newly formed and returned by star formation; the 
second term represents the metal mass locked into long- 
lived stars; the third term represents the metal mass in 
the ISM that is carried away by outflows. 

We note that the chemical yield y and the return frac¬ 
tion R are approximately constant parameters, which 
are largely determined by the stellar initial mass func- 
tion (IMF) and weakly d epend on metallicity ( e.g. BC03 
Bruzual fc Charlotll2Ci03ll . For a Chabrier IMF l)Chabrierl 
2003(1 . R = 0.46 is relevant to adopt for the instantaneous 
recycling approximation. The outflow mass-loading fac¬ 
tor ?7 is a variable to be constrained in this paper. In 
principle, the outflow mass-loading factor 77 is a variable 
that can change with time and from galaxy to galaxy. In 
this paper, we assume that 77 can arbitrarily vary as a 
function of the present-day galaxy stellar mass, but is a 
constant for a given galaxy (i.e. is not a function of time). 
This is equivalent to assuming that the parameter 77 is 
an effective mass-loading factor averaged over the history 
of the galaxy and weighted by the star formation rate, 
i.e. 77 = f*° y(t)(/>(t)dt/ c/)(t)dt, where is the time 

dependent star formation rate. We note that because 
stellar mass can only increase for an isolated galaxy, one 
can use stellar mass as a clock to integrate the increase 
of the metal mass over time. 

Using this simple model, we can derive a set of equa¬ 
tions to describe the relationship between metallicities 
and parameters characterizing star formation and out¬ 
flow. In the first set of derivations, we seek the rela¬ 
tionship between the stellar-phase metallicity, gas-phase 
metallicity and the outflow mass-loading factor. The 
change of the metal mass locked into long-lived stars in a 
time interval is the stellar mass formed in the time inter¬ 
val multiplied by the instantaneous gas-phase metallicity 
Zg, as 

dMz_* = ZgdM* . (2) 

From this equation, one realizes that the term ZgdM* in 
Equation o is just the metal mass that is locked into 
long-lived stars. The integral of Equation (HD immedi¬ 
ately yields that 


Mz_g = 7/M* 



1-R 


Z^cM^, , 


( 3 ) 


where Z* is the averaged stellar-phase metallicity. 
Adopting the definition for the gas-phase metallicity, 
Zg = Mz_g/Mg, we hnd that the gas-phase metallicity 















































































3 


is 


Z, = 


y 



l-R 




( 4 ) 


where rg is the gas-to-stellar mass ratio of galaxies at the 
present day, e.g. rg = Mg/M*. We note that how the 
gas-to-stellar mass ratio evolves in the past is irrelevant 
here, because this equation only describes the present- 
day metallicities of a galaxy. This equation shows that 
the gas-phase metallicity decreases with an increasing 
gas-to-stellar mass ratio, as the cold gas dilutes the 
metallicity of the ISM. It also shows that the gas-phase 
metallicity decreases with an increasing outflow mass¬ 
loading factor, as outflow expels metal-enriched gas out 
of galaxies. We can rewrite the equation as follows to 
determine the mass-loading factor using metallicities. 


r; = (1 - i?) 


JL 



( 5 ) 


This equation shows that, for a given galaxy, higher 
metallicities and higher cold gas-stellar mass ratio indi¬ 
cate a lower outflow mass-loading factor. This is easy to 
understand because strong outflow will reduce the metal 
content and gas content of the galaxy. 

In the second set of derivations, we adopt the gas-phase 
metallicity at different redshifts as data constraints. To 
do so, we need to assume how the cold gas metallicity 
evolves with stellar mass growth. Here we first adopt a 
simple assumption that the gas-phase metallicity follows 
a power-law trajectory in the Zg-M* plane as the galaxy 
evolves, e.g. 




fM4t) 

'v M* 


A* 


( 6 ) 


in which the power index, is assumed to vary as a func¬ 
tion of galaxy stellar mass, and is to be determined. The 
variation of /i reflects different ways in which galaxies can 
evolve their metallicities as their stellar masses increase. 
For example, inflow can decrease metallicity, and out¬ 
flow can increase or decrease the metallicity depending 
on how metals are loaded in the wind. The assumption 
of power-law trajectories allows galaxies with different 
masses to evolve along different paths. Based on this 
simple assumption, by integrating Eq[T]we find that 


M,,g = yM* - (^1 + • ( 7 ) 

Applying the definition of the gas-phase metallicity, we 
can rewrite the equation to be an expression for the mass¬ 
loading factor, 


11 = (l-R) 





( 8 ) 


This equation shows that we can also constrain the av¬ 
eraged mass-loading factor 77 , when the trajectory of a 
galaxy in the Zg — M* plane is determined. 

In the above derivations, we have assumed that all the 
newly produced metals are returned into the ISM, so that 
the yield affecting the chemical evolution of a galaxy 
equals to the intrinsic nucleo-synthesis yield. Because 
SN ejecta can preferentially transport over-enriched (rel¬ 
ative to the mean ISM) gas out of the potential well of 


a galaxy as shown i n hydro-dynamical simulations (e.g. 
iCreasev et al.l[2015t iMelioli et al.ll2015fl . it is reasonable 
to assume that the “retained” yield, y, is less than the 
intrinsic yield, y. In the following analysis, we vary the 
retained yield to explore this effect. 


3. DATA 


The derivations presented in the last section demon¬ 
strate that one can draw inferences about star forma¬ 
tion and feedback when certain observational data are 
given. Based on our derivations, these observational 
data include the stellar-phase metallicity as a function 
of galaxy stellar mass, the gas-phase metallicity as a 
function of galaxy stellar mass of local galaxies, and of 
galaxies at different redshifts. To demonstrate the con¬ 
straining power of these observational relations based 
on the analytic derivations shown in the previous sec¬ 
tion, we choose to use the mean relation of each of ob¬ 
servational results. For the stellar -phase metallic i ty, we 
ado pt the observationa l results of IGallazzi et al.l (120051 ) 
and iKirbv et al.l (I2013D to cover a wide range of stellar 


masses between ^ lO^M^ and lO^^M^. As shown by 
IKirbv et al.l (I2013I1 . these two pieces of data join with 
each other remarkably well if the relation determined by 
the Kirby data is extrapolated to high stellar masses. 
In the analysis of this paper, we adopt a function that 
follows the two data sets and smoothly joins them as 
shown in Fig. [TJ We note that because the uncer¬ 
tainty in the stellar metall icities at the low mass end 
of the IGallazzi et al.l (120051) relation is large, the small 
deviation between their mean relation and our adopted 
function is ignored in this paper. In addition, there are 
two biases in the stellar-phase metallicity. First, the ob¬ 
served metallicity is weighted by young and luminous 
stars, which tend to be metal rich. Second, the stellar- 
phase metallicity measurements are mainly sensitive to 
iron abundances rather than the a elements, which are 
the species measured in the gas-phase m etallicity mea¬ 
surements. Following iPeeoles et al.l (I20I4D . we have cor¬ 
rected these biases using their Eq.(7) and Eq.(AI), taking 
into account the different solar metallicity adopted in the 
reference paper. Eor the gas-phase metallicity relations, 
there are a larger number of results published in the liter¬ 
ature. While the observed gas-phase metallicity relations 
have small scatter at a given redshift (typically 0.1 dex), 
they have substantial differences in their normalizations 
and shapes owing to large uncertaint ies in the calibra¬ 
tion of the metallicity measurements (jKewlev fc EllisonI 
120081) . To take these uncertainties into account in our in¬ 
ferences, we adopt multiple observational results. Eor the 
gas-phase metallicity-stellar mass relati on of local star 
formin g galaxies, we ad o pt th e results ofiTremonti et al.l 
(120041): iM aiolino et al.l (120081) : iZahid et al.l ( 20131) and 
lAndrews fc Martinil (|2013l) for galaxies with stellar mass 
in a range between ~ 10 ®Mq and lO^^M^, and lLee et al.l 
(|2006l) for lower mass galaxies in the range between 
^ lO^Mg and lO^Mg. Eor higher redshifts, we adopt 
th e metallicity-ste l lar m ass r elation at z ^ 2.2 compiled 
bv IMaiolino et al.l (I2008D and IZahid et al.l (120131) , respec¬ 
tively. In additio n, we also us e the ga s-phase metallicity 
measurements of lHenrv et al.l (j2013al ftJ) for a lower mass 
range between 10 ®Mq and lO^^M^ as a complement. To 
infer the mass-loading factor with the approaches demon- 
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strated in the previous section, one also needs to know 
the gas-to-stellar mass ratio of galaxies. We adopt a fit¬ 
ting formula for the mean relation o f star forming galax¬ 
ies compiled bv lPeenles et al.l l)2014[ l. namely 

logrg = -0.48 log +4.39. (9) 

As shown by the authors, this fitting formula captures 
the mean relation between the gas mass ratio and galaxy 
stellar mass well and agrees very well with many observa¬ 
tiona l estimates (e.g. lMcGaughl[200M 120121 : iLerov et’aTI 
120081 : iPaoastergis et al.l 1201211 . Unlike the metallicity- 
stellar mass relations, the cold-gas to stellar mass ra¬ 
tio has a significantly large scatter, ~ 0.5dex. We add 
this scatter into the gas mass ratio as upper and lower 
bounds with A log rg = 0.5 to demonstrate how this scat¬ 
ter propagates into our results. For all the dat a sets, we 
have c orrected the stellar mass by assuming a iChabri^ 
(1^)1 IMF. 



Fig. 1.— Observational results of metallicity as a function of 
galax y stellar mass . The stell ar-phase met a llicity relations are 
from iGallazzi et al ] 1)20051') and IKirbv et al.l ((2^Q), denoted by 
the solid lines of G05 and K13, respectively. The green dash- 
dotted line shows a smooth function joining both of the data 
sets. The loca l galaxy gas-phase metallicit y relations a r e from 
Tremonti et al.l 1120041b iMaiolin o et al.l 1120081). IZahid et al.l 1120131b 
^idrew^^Martinil ^20131) . and ILee et al.l feood) . which are de¬ 
noted by colored solid lines of T04, M08, Z13, A13, and LOG, re- 
s pectively. The gas-pha se metallicity relations at higher redshifts 
of IMaiolino et all 11200 80 and IZahid et al.l II2013I) are shown by long- 
dashed lines for 2 : 0.7 and sho rt-dashed lin e s for 2 : 2.2. The 

gas- phase metall i city re lations of lHenrv et al.l l)2013al ') for 2 0.6 

and fHenrv et all Il2013bl) for 2 : 1.8, denoted by H13, are shown by 

long-dashed and short-dashed orange lines. The horizontal short 
dashed line denotes the intrinsic yield y assumed in the paper. 


4. RESULTS 

4.1. Outflow mass-loading factor Inferred from local 
metallicity relations 

As we demonstrated in ^ without assuming any par¬ 
ticular star formation history or metallicity enrichment 
history, our simple model can match the stellar-phase 
and gas-phase metallicity relations at a given redshift and 



Fig. 2.— The constrained outflow mass-loading factor 77 derived 
from different combinations of observational data. The shaded 
regions and the solid lines are the results inferred from the stellar- 
phase metallicity and the gas-phase metallicity of local galaxies. 
The shaded regions show rj obtained by assuming that all pro¬ 
duced metal mass is mixed into the ISM but varying the cold gas 
to stellar mass ratio. The upper and lower bounds encompass the 
variations resulted from the ±0.5dex scatter of the cold gas mass 
ratio for given stellar mass. The solid lines show the results assum¬ 
ing half of the produced metal is directly expelled from galaxies, 
i.e. y = 0.5y with the fiducial cold gas mass ratio. The dashed 
lines denote the results inferred from the gas-phase metallicity re¬ 
lation at multiple redshifts (see ED- Different colors denote dif- 
ferent data sources. The brown line represents the result using 
lAndrews Sz. Martinil II2013I) data. The purple lines denote the re- 
sults using data. The blue line shows the 

result using ITremon^ et al.l 1I2004I~) data. The cyan lines show the 
results using IZahi^^^aLI (120131) data. The orange line shows the 
results usin gfHe nrv^ aj4(l2Q13alfblj data. The green line shows the 
result using [T7ee et al.l 1)20061') data. The black dotted and dashed 
lines denote the upper limit estimates of the mass-loading factor 
allowed by two different re-ionization models. 


can constrain the mass-loading factor. In this subsec¬ 
tion, we adopt the combined the stellar-phase metallicity 
as a function of stellar mass and the g as-phase metal¬ 
licity r esults of iTremonti et alJ (I20Q4D: iMaiolino et alJ 
(1200811: IZahid et ah ( 2013D : lAndrews &: Martinil ( 2013D . 
and ILee et al.l (|20 'q 5) of local galaxies to infer the out¬ 
flow mass-loading factor, 77, as a function of galaxy stel¬ 
lar mass. We also inc l ude th e gas-phase metallicity re¬ 
lation of iHenrv et al.l (j2013a[) at z = 0.6 in this analy¬ 
sis. We choose each of the gas-phase metallicity relation 
and combine it with the stellar-phase metallicity rela¬ 
tion and the gas-to-stellar mass ratio to feed into Equa¬ 
tion ©■ At first, we assume all of the metals produced 
are retained in the galaxy and mixed into ISM before 
they can be expelled by outflows, i.e y = y. In other 
words, outflows have the same metallicity as the ISM 
with no further metallicity enhancement. We vary the 
normalization of the gas-to-stellar mass ratio, rg, up and 
down by 0.5 dex to encompass the random variation of 
the gas mass for given stellar mass. The resulting mass¬ 
loading factor as a function of galaxy stellar mass for 
each data set is shown by the shaded bands in Fig. [51 
For the all results, the general trend is that the mass¬ 
loading factor is higher for low-mass galaxies and de¬ 
creases rapidly for higher stellar masses. As one can see. 
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the inferred r]-M^ relations by using different gas-phase 
metallicity relations are very similar to each other at the 
high-mass end, where the inferred mass-loading factor 
is less than 2. Another interesting feature is that the 
inferred mass-loading factor tends to flatten out at the 
very high mass end. The reason for this is that mergers 
become increasingly important for higher mass galaxies, 
and the mass-loading factor we constrain is considered to 
be an average over all progenitors, which have relatively 
lower masses. In the lower mass regime, the scatter in 
the gas mass ratio produces increasingly larger variations 
in the inferred outflow mass-loading factor. Neverthe¬ 
less, the upper bound for ry is well defined. For galaxies 
M* = lO^Mg, the upper-bound mass-loading factor in¬ 
ferred from all the data combinations is less than 20. 
For the same stellar mass, galaxies with higher gas frac¬ 
tion may have significantly lower rj. In addition, different 
gas-phase metallicity results have different level of sen¬ 
sitivity to the variation of the gas mass ratio. A higher 
gas-phase metallicity not only suggests a systematically 
lower mass-loading factor, but also a more sensitive de¬ 
pendence on the gas mass ratio. Second, we assume that 
SN ejecta carry away half of the produced metal mass 
without mixing it into ISM. In this situation, y = 0.5t/, 
and we show the inferred rj with the median gas mass 
ratio by solid lines in Fig 121 When a fraction of metal 
leaks from galaxies without mixing into ISM, the ob¬ 
served metallicity relations yield a much lower rj. At the 
high-mass end, all the data combinations suggest that 
the mass-loading factor is below 1. At M* = lO^M^, the 
mass-loading factor is below 10, and it increases with de¬ 
creasing stellar mass. At the low-mass end of the mass 
range lO’^Mg), the mass-loading factor is about 20. 
It is worth noting that we chose to use a rather high 
metal yield, y = 0.07, which is about five times of the 
Solar metallicity, and at the high end of usual values pre¬ 
dicted by neucleo-synthesis models (|Wooslev fc Weaven 
119951: iNomoto et alJ[T^97bllall . We chose this high yield 
parameter to demonstrate an upper limit for the mass¬ 
loading factor. If a lower metal yield is chosen instead, 
the inferred mass-loading factor decreases. 

For 1 0^M,„ gala xies, the mass-loading factor inferred 
bv the iL ee et all J2006f ) gas-phase metallicity and the 
iKirbv et al.l (1201^ stellar-phase metallicity relations can 
be as high as ^ 20. We note, however, that low-mass 
galaxies are expected to be hosted by low-mass halos, 
and re-ionization would prevent a large fraction of bary- 
onic mass from collapsing into such low-mass halos in 
the first place, so that galaxies forming in those halos 
would never have contained as much baryonic mass as 
these large mass-loading factors implied, unless the out¬ 
flow materials are rapidly reaccreted back into the galaxy. 
We note that ry is the averaged net mass-loading factor, 
while the instantaneous mass loading in individual star- 
burst galaxies may be higher if there is significant baryon 
mass recycled from early times. Under the assumption 
that reaccretion is not important, the total baryonic mass 
in stars, cold gas, and ejected by outflow should be lower 
than the total baryonic mass that can collapse into the 
host halo, i.e. 


M* -I- M*rg -I- ?y 


I -i? 


A /b/reion(Afvir5 ^)Afvir i 


( 10 ) 


where /b = 0.17 is the cosmic baryon fraction, and /reion 
is the fraction of baryon mass that can collapse into 
a halo with virial mass Mvir at redshift ^ due to re¬ 
ionization. This inequality yields that 


y<{l-R) 


/b/n 


Mv 




- re - 1 


( 11 ) 


To estimate the limit for the mass-load i ng fa ctor, we 
adopt a re-ionization mod el bv iGnedinI (1200011 with a 
fitting formula proposed bviKravtsov et ahT i 20041 1 and a 
recent model bv INoh &: McQuinnI (|2014ir to compute the 
total baryon mass fraction as a function of halo mass at 
z = 0. We then adopt the abundance matching model of 
iBehroozi et al.l (l2013f) and extrapolate it to the low-mass 
end (M* ^ lO^Mg) to determine the stellar mass-halo 


mass relation. Using these and Equation dm , we can 
derive an upper limit for the mass-loading factor allowed 
by re-ionization. The thick dashed line and dotted line in 
Fig [2] show the upper limits set by the Gnedin model and 
the Noh & McQuinn model, respectively. In both cases, 
the limiting outflow mass-loading factor drops sharply 
when galaxy stellar mass goes below ~ 10 ®Mq. Be¬ 
low this mass scale, re-ionization prevents baryons from 
collapsing into the low-mass halos hosting these galax¬ 
ies. These low-mass galaxies should therefore never have 
contained as many baryons as is implied by their inferred 
mass-loading factors given their existing stellar mass and 
cold gas mass. To reconcile the observed metallicities and 
cold baryon masses of low-mass galaxies, a more plausible 
model seems to require metal-enriched outflows, which 
carry away a large fraction of metals directly from SN 
ejecta without mixing much metals win ISM. In this sce¬ 
nario, the retained metal yield can be lower than we used 
here, yielding lower mass-loading factors. 


4.2. Outflow mass-loading factor inferred from 
gas-phase metallicity relations at different redshifts 

We now use the gas-phase metallicity relation at mul¬ 
tiple redshifts to carry out the inference described in Sj2l 
We adopt th e gas-phase met a llicity -stellar mass relation 
compiled by iMaiolino et al.l (120081) for galaxies at two 
different re dshifts z = 0.07 an d z = 2.2, the relation 
compiled bv iZahid et al.l (l2013fl at )g = 0.08 and z = 2.3, 
and the relation measured bv iHenrv et all (j20I3alfO ) at 
z = 0.6 — 0.7 and z = 1.3 — 2.3 for lower stellar masses 
(A/, = 10® — 10 ^°Mq ). As we described in we use this 
data to determine the trajectories along which galaxies 
evolve in the Z-M^ diagram. For the simple case we 
demonstrated in 321 we essentially need to determine the 
logarithmic slope, y, as a function of galaxy stellar mass 
using the observed gas-phase metallicity relations at two 
separate redshifts. One way to determine /r(M*) is to 
use a realistic star formation history of a galaxy with 
a stellar mass A7, at a lower redshift to determine the 
stellar mass of its typical progenitor at a higher redshift, 
and then determine where the progenitor galaxy is in the 
Zg-M^ diagram at the higher redshift by using the gas- 
phase metallicity-stellar mass ration at that redshift. We 
adopt a fitting model for the star formati o n hist ories of 
star forming galaxies proposed bv iLeitn^ ( 120121 ). 

{1 + zr, ( 12 ) 
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Fig. 3.— The grey lines show three examples of the power-law 
trajectories of galaxies in the Zg-M* diagram. These trajectories 
are determined by ass uming star formation histories parameter¬ 
ized by Equation GUi an d the gas-phase meta llicity relations at 
= 0.07 and 2 : = 2.2 from lMaiolino et al.l II2008I) . The blue dashed 
line shows the predicted gas-phase metallicity as a function of stel¬ 
lar mass at 2 = 0.7 using the deter mined trajector i es. T he solid 
lines are the observational results of IMaiolino et aT1 1120081) at dif¬ 
ferent redshifts as noted in the legend. The inserted diagram shows 
the inferred /i (the power index for the metallicity-stellar mass tra¬ 
jectory) as a function of final stellar mass in units of M^. 


where parameters Aq = 3.24MQyr“^, a = 3.45, and 
P = —0.35 are determined by matching the abundances 
of star forming galaxies out to high redshift (z ^ 4). Us¬ 
ing this model, we can integrate the star formation rate 
history to predict the stellar mass at any given redshift 
for a galaxy with a given final stellar mass. For simplic¬ 
ity, we again assume instantaneous recycling when inte¬ 
grating SFR histories. For more accurate models taking 
into acc ount time-dependent recyc li ng, readers ar e re- 
ferred to lLeitner &: Kravtso^ (1201111 . ILeitned (|2012fl and 
iLu et al.l ( 2015all . Assuming a galaxy evolves along a 
power-law trajectory on the gas-phase metallicity-stellar 
mass plain between two different redshifts, we can find 
where the trajectory determined by the star formation 
history intersects with the gas-phase metallicity-stellar 
mass relations at the different redshifts. Using the 
two intersection points in the diagram, we 

can determine the logarithmic slope as /r = (log 2° — 
logZg)/(logM° - logM^^), where the superscripts “0” 


and “1” denote quantities at two different redshifts of a 
galaxy. FigI3] shows an example of how the /r param¬ 
eter is_detemnri^ as a function of stellar mass using 
the IMaiolino et all (1200811 data. In the figure, we show 
three trajectories with different final galaxy masses in 
grey lines and o ver-plot them w i th the gas-phase metal¬ 
licity relation of IMaiolino et al.l (|2008ll at z = 0.07, 0.7, 
and 2.2. For each given final stellar mass at z = 0.07, we 
compute the stellar mass at z = 2.2 using the star for¬ 
mation history of EqlT^l We read the gas-phase metal¬ 
licity of the galaxy with the stellar mass from the gas- 
phase metallicity relation at the corresponding redshift. 
The open circles along a grey line are the positions of 
the galaxy at the two redshifts, z = 0.07 and 2.2. The 


grey lines connecting the two circles are the power-law 
trajectories of the three example galaxy masses. The 
grey squares between the two circles on each trajectory 
mark the predicted gas-phase metallicity-stellar mass re¬ 
lation at z = 0.7. We predict those quantities at the 
redshift by interpolating along the star formation his¬ 
tory and the power-law Z~M^ trajectory. Connecting 
these predictions, the dashed line shows the predicted 
gas-phase metallicity-stellar mass relation at z = 0.7. 
As one can see, the interpreted gas-phase metallicity re¬ 
lation at z = 0.7 agrees with the observational result at 
the same redshift remarkably well, which demonstrates 
that the power-law trajectory sufficiently well captures 
the evolution of a galaxy in the Zg-M* plane. In the 
inserted panel in FiglSl we show the determined ^ as a 
function of final galaxy stellar mass. The resulting ^ 
slowly decreases from 0.73 at M* = lO^M^ to 0.51 at 
M, = lO^^Mg. IPeenles fc Somervillel (I20li^ have per¬ 
formed more detailed inference on the gas-phase metallic¬ 
ity trajectories by assumin g a fundamental meta llicity- 
stellar mass-SFR relation (|Mannucci et al.l 1201011 with¬ 
out enforcing the trajectories to follow a power -law form. 
From Figure 1 of iPeeples fc Somervillel (|2013f) . one can 
find that their trajectories are nearly power-law over a 
large range of redshifts, which in turn supports our sim¬ 
ple assumption of power-law trajectories made in Equa¬ 
tion ([6]). 

Combining the determined /r(M*) and the data and 
inserting these terms into Equation ([5|), we can again 
compute the mass-loading factor as a function of stel¬ 
lar mass. The results are shown as long dashed lines 
in Fig[^ As one can see, the results obtained from 
this approach are qualitatively similar to those inferred 
without assuming a star formation history in M.ll The 
constrained mass-loading factor decreases with increas¬ 
ing stellar mass. Quantitatively, using the same yield, 
y = y = 0.07, the mass-loading loading factors inferred 
in this approach are a factor of 2-3 lower than those in¬ 
ferred in (14.11 The difference may be due to the calibra¬ 
tion of different data sets or inconsistency between the 
empirical star formation history and the real star forma¬ 
tion history of the star-forming galaxies included in the 
gas-phase metallicity measurements. 

4.3. Metallicity history 

The metallicity of a star approximately tells us the 
metallicity of the ISM at an early time when the star 
formed, while the gas-phase metallicity tells us the cur¬ 
rent metal content of the ISM. Using these two pieces 
of information, one can gain some insight into the metal 
enrichment history of galaxies. The left panel of Fig¬ 
ure in shows the ratio between the stellar-phase metal¬ 
licity and the gas-phase metallicity as a function of stel¬ 
lar mass of local galaxies. For the gas-phase metallic¬ 
ity rel at ions, we adopt the resu lt s oflAndrews &: Martinil 
(I20I3I1. IMaiolino et al.l (1200811. iTremonti et al.l (1200411 . 
IZahid et al.1 ( 20I3I1 . and lLee et al.l (|2006ll for local galax- 
ies. For th e stell ar-phase m etallicity, w e use the 
iGallazzi et al.l (|2005l l result and iKirbv et al.l (|20I3[1 re¬ 
sult for relevant stellar mass range. Although these data 
sets have different amplitudes due to different calibra¬ 
tions, they show a systematic trend that low-mass galax¬ 
ies (M* < 10 ^°Mq) have a relatively lower stellar-phase 
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Fig. 4.— The left panel shows the stellar-phase metallicity to gas-pha s e metallicity rat i o as a f unction of st e llar m a ss. The local galaxy 
gas-ph ase metallicity relatio ns are adopted from ITremonti et al] Il2004l 'l , IMaiolino et al.l <1200811 , IZahid et al.l Il201,'ll l , lAndrews fc Martinil 
1 I 2 OI 3 I I. and lLee et al.! <120061 1. which are denoted by colored solid lines of T04, M08, Z13, A13, and L06, respectively. Either the Gallazzi et 
al. result (G05) or the Kirby et al. result (K13) for the stellar-phase metallicity is adopted for the relevant stellar mass ran ge. The right 
panel shows the stellar-phase metallicity to gas-phase metallicity ratio as a function of parameter {i and fc in Equation insii. 


metallicity to gas-phase metallicity ratio than the high- 
mass ones. This trend indicates that galaxies with dif¬ 
ferent masses follow different ways of metal enrichment. 

Assuming that the gas-phase metallicity trajectory in 
the Zg — M* plane follows a power-law function, we find 
that the ratio of the stellar phase metallicity and the 
gas-phase metallicity is 


trajectory of Zg when M* > Me, and Me is a character¬ 
istic mass scale at which the Zg trajectory switches from 
an exponential function to a power-law. Below Me, the 
gas-phase metallicity increases exponentially with mass. 
This model captures the idea that feedback might be 
powerful enough to expel all produced metals at the early 
epoch of galaxy formation, and it becomes less powerful 
at late time when the galaxy mass is high so that metals 
can be retained. 

Using this model for the trajectory, we find that 


where fj, is the power-low index of the trajectory. What 
this equation tells us is that a lower stellar-phase metal¬ 
licity to gas-phase metallicity ratio corresponds to a 
steeper gas-phase metallicity trajectory (larger fi). The 
left panel of Figure H] shows that high-mass galaxies have 
higher Z*/Zg, which means they need to have lower 
/r than low-mass galaxies. This is consistent with the 
/r — M* relation we derived from the gas-phase metallic¬ 
ity relation at different redshifts in H4.2I (see the inserted 
panel of Fig.[3l). This suggests that high-mass galaxies 
seem to maintain a slow increase of metallicity for a long 
time in the past, while low-mass galaxies have rapidly 
enriched their metallicities in recent times. 

We can also use a more sophisticated model to describe 
the trajectory, as the power-law model might be too sim¬ 
ple to capture the entire evolution, especially at early 
times. We choose the following model using a power-law 
to describe the late time evolution, and an exponential 
term to capture the possible rapid increase of metallicity 
following a starburst phase at the early time, i.e. 


ZS) 


Zg,o 

K 



exp 


Me 

M* (t) 


(14) 


where K = exp is a normalization fac¬ 

tor, M*.o and Zg^Q are the final stellar mass and gas-phase 
metallicity, /r parameter is the power-law slope of the 


^ =/cexp(/c)r(-^-1,/c), (15) 

where fe = and r(s,a:) is the upper incomplete 

gamma function: F(s,a;) = We show the 

two-dimensional function in the right panel of Figure |4l 
The contours in the figure show the combinations of fe 
and fj. that give rise to a constant Z*/Zg. As one can 
see, when Z*/Zg ^ 1, the gas-phase metallicity track is 
required to be flat (^ ~ 0 and fe is very small). When 
Z^jZg becomes smaller, the gas-phase metallicity trajec¬ 
tory has to be a steep function of stellar mass, either with 
a large power-law slope /i or a higher fe. To result in this 
type of metallicity trajectory, a large fraction of stellar 
mass needs to form in an early starburst with very low 
metallicity, and to retain metals in the ISM to rapidly 
enrich the gas-phase metallicity at late time. 

This inference has two implications. First, the low 
ZtfjZg ratio of low-mass galaxies seems to indicate that 
a relatively large fraction of stellar mass of these galax¬ 
ies formed in early starbursts when the ISM metallic¬ 
ity is still low. Second, low-mass galaxies, as opposed 
to high-mass galaxies, are more rapidly enriching their 
metallicity, indicating a downsizing effect that low-mass 
galaxies are less evolved than the high-mass ones also in 
metallicity. 


5. DISCUSSION 






































In this paper, we have introduced a simple analytic 
model to follow the metallicities of galaxies. When the 
model is matched to the gas-phase and stellar-phase 
metallicities as functions of stellar mass of galaxies at 
different redshifts, these observational data in turn can 
constrain star formation and feedback. 

With minimum assumptions about how metals pro¬ 
duced by star formation are mixed into the ISM and 
how outflow can affect the baryonic content of galaxies, 
we demonstrate that metallicities of galaxies provide use¬ 
ful constraints on the strength of galaxy outflows. Using 
the model, we have derived two approaches to constrain 
the outflow mass-loading factor as a function of galaxy 
stellar mass. In the first approach, we find the stellar- 
phase metallicity, gas-phase metallicity and the cold-gas- 
to-stellar mass ratio as functions of galaxy stellar mass 
at a given redshift can already provide strong constraints 
on the upper limit for the mass-loading factor. The 
fundamental reason that the combination of these data 
can constrain outflow is because the gas-phase metallic¬ 
ity represents the metal content of galaxies at a fixed 
time, and the stellar-phase metallicity represent a time- 
averaged value for the metal content of galaxies. The 
combination of these two quantities provide constraints 
on how the metal content changes over time when galax¬ 
ies build up their stellar mass. The metal mass that is not 
found in galaxies should have been carried away by SN 
ejecta or galactic outflow. The second approach relies on 
the gas-phase metallicity relation at different redshifts. 
Once the manner in which galaxies change their metal¬ 
licity as they increase their stellar mass is determined, 
the metallicity relations at different epochs can also con¬ 
strain the strength of outflow. 

Combining multiple observational measurements for 
the stellar mass-metallicity relations, we have derived the 
outflow mass-loading factor based on the simple model. 
The results show that, in spite of the large uncertain¬ 
ties in the metallicity measurements, the inferred out¬ 
flow mass-loading factors have a well defined upper limit. 
Even with generous (high) metal yield and assuming all 
produced metals are mixed into the ISM, the inferred 
mass-loading factor cannot be higher than 20 for lO^Mg 
galaxies and it drops quickly for higher galaxy masses. 
On the very low mass end, where re-ionization is effec¬ 
tive to block a significant fraction of baryons from col¬ 
lapsing into the galaxies, strong outflow is not a plau¬ 
sible explanation for the cold baryon mass fraction and 
the metallicity relations. To explain the data, we find 
that a significant fraction of produced metals need to 
be ejected directly out of the galaxies with SN ejecta 
without mixi ng in the I SM, w hich is consistent with 
the results of iDalcantonI (|2007ll . who found that metal- 
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enriched outflow is the only viable mechanism to repro¬ 
duce the low retained yields of gas-rich low-mass galax¬ 
ies. Observations of metal content of winds from dwarf 
starburst also suggest that almost all the metals pro¬ 
duced i n starbursts a.re eje cted directly out of the galax¬ 
ies le.g. lMartin et al.ll2002[ l. providing observational sup¬ 
port of metal-enriched ejection. If a moderate fraction 
(half) of metal mass is involved in the metal-enriched 
ejection, the mass-loading factor is expected to be at 
most only a few (77 < 8 ). This low mass-loading fac¬ 
tor conflicts with many galaxy formation models, which 
require much higher mass-loading factors, especially for 
low-mass galaxies (M* < I 0 ^°M ^), to exp lain low-mass 
end o f the stellar mass function (jLu et ani2ni4 iBensonI 
[Ml . In many successful galaxy formation models, 
strong outflow is needed because baryons other than ob¬ 
served stars and cold gas cannot stay in the halos with¬ 
out forming stars. The only way to keep low-mass ha¬ 
los from having a too high baryon mass fraction is to 
eject a large amount of baryonic mass out of the halo. 
What the metallicity relations seem to show in this pa¬ 
per, however, is that when galaxies increase their stel¬ 
lar mass (and simultaneously produce metals), the metal 
mass in galaxies increases substantially. The increase in 
the metal mass suggests that there should not be very 
strong outflow to take away metals. 

The data show that low-mass galaxies tend to have a 
lower stellar-phase to gas-phase metallicity ratio. Using 
simple analytic models, we conclude that this trend sug¬ 
gests that low-mass galaxies tend to increase their gas- 
phase metallicity rapidly in late times, while high-mass 
galaxies have gone through the rapid metal enrichment 
phase in early times, which indicates an imprint of the 
downsizing effect in the metallicity evolution of galax¬ 
ies. The analysis also suggests that before the rapid en¬ 
richment epoch, a fraction of stellar mass is formed at 
early times when the gas-phase metallicity was low. This 
is consistent with recent observational results that local 
dwarf galaxies have more than hal f of their stellar m ass 
form in early phase prior to z = 2 (|Weisz et al.ll2014ll . 

Finally, we stress the importance of the accuracy of 
the metallicity relations and cold gas mass fraction in ob¬ 
servational measurements. We have adopted simplified 
version of existing observational results to demonstrate 
the constraining power of the data. Better observational 
determination of these scaling relations will definitely im¬ 
prove the constraints when more sophisticated model in¬ 
feren ces. such as Bayesian model inference le.g. lLu et al.l 
unni), are adopted. 
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